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We present the next-to-leading-order (NLO) QCD corrections to the inclusive total cross section for the pro¬ 
duction of a Higgs boson in association with a top anti-top quark pair within the Standard Model at the Tevatron 
and the LHC. 


1. Introduction 

The search for the Higgs boson of the Stan¬ 
dard Model (SM) is one of the major tasks of 
the next generation of high-energy collider ex¬ 
periments. The direct limit on the SM Higgs 
boson mass, Mh, from LEP2 searches and the 
indirect limit from electroweak precision data 
strongly suggest the existence of a light Higgs 
boson, Mh < 193 GeV (95 % C.L.) [1], which 
may be within the reach of the Fermilab Teva¬ 
tron pp collider [2]. However, since the domi¬ 
nant SM Higgs production channels are plagued 
with low event rates and large backgrounds, the 
Higgs boson search at the Tevatron requires high 
luminosity, and all possible production channels 
should be considered. The production of a SM 
Higgs boson in association with a top-quark pair, 
pp —> tth , can play a role for almost the entire 
Tevatron discovery range [2,3]. Although the tth 
event rate is small, the signature of such events is 
quite spectacular (W + W~bbbb). 

At the CERN LHC pp collider, the associated 
tth production mode will play a crucial role in 


the 110 GeV< Mf, <130 GeV mass region both for 
discovery and for precision measurements of the 
Higgs boson couplings. This process will provide 
a direct measurement of the top-quark Yukawa 
coupling and will be instrumental in determin¬ 
ing ratios of Higgs boson couplings in a model 
independent way [4]. Such measurements could 
help to distinguish a SM Higgs boson from more 
complex Higgs sectors, e.g., as predicted by su¬ 
persymmetry, and shed light on the details of the 
generation of fermion masses. 

As for any other hadronic cross section, the 
next-to-leading-order (NLO) QCD corrections 
are expected to be numerically important and are 
crucial in reducing the (arbitrary) dependence of 
the cross sections on the factorization and renor¬ 
malization scales. Here we summarize the cal¬ 
culation of the NLO QCD predictions for associ¬ 
ated tth production at the Tevatron and the LHC. 
Results for the Tevatron have been presented in 
Refs. [5-8] and for the LHC in Refs. [7-9]. The 
results of the two groups are in very good agree¬ 
ment within the statistical errors. 
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2. QCD corrections to tth production 

The inclusive total cross section for pp (or pp) 
—> tth at 0(ag) can be written as: 

o- N lo(pp —> tth ) = 1 [ dxidx 2 (1) 

U 1 + S V ) 

■[J 7 ? (X!, n)ff p) (x 2 ,p)d” LO {ii) + (l<-> 2)] , 

where T p ' v are the NLO parton distribution func¬ 
tions for parton i in a proton/antiproton, defined 
at a generic factorization scale fif=p,, and a^ LO 
is the 0(a%) parton-level total cross section for 
incoming partons i and j, renormalized at an ar¬ 
bitrary scale fi r . We take p r = p,f = p. The 
NLO parton-level total cross section ij 5? [0 con¬ 
sists of the 0(a 2 ) Born cross section and the 
0(a s ) corrections to the Born cross section, in¬ 
cluding the effects of mass factorization. It con¬ 
tains virtual and real corrections to the parton- 
level tth production processes, qq —» tth and 
gg — > tth, and the tree-level ( q , q)g initiated pro¬ 
cesses, ( q,q)g —> tth(q,q), which are of the same 
order in a s . The main challenges in the calcu¬ 
lation come from the presence of pentagon dia¬ 
grams in the virtual corrections with several mas¬ 
sive external and internal particles, and from the 
computation of the real part in the presence of 
infrared singularities. 

2.1. Virtual Corrections 

The calculation of the virtual corrections to the 
qq initiated sub-process is described in detail in 
Ref. [6]. The calculation of the virtual correc¬ 
tions to gg —> tth is technically similar. The ba¬ 
sic method is to reduce each virtual diagram to a 
sum of scalar integrals, which may contain both 
ultraviolet (UV) and infrared (IR) divergences. 
The finite scalar integrals are evaluated by us¬ 
ing the method described in Ref. [10] and cross 
checked with the numerical package FF [11], The 
scalar integrals that exhibit UV and/or IR diver¬ 
gences are calculated analytically. Both the UV 
and IR divergences are extracted by using dimen¬ 
sional regularization in d= 4—2e dimensions. The 
UV divergences are then removed by introducing 
a suitable set of counterterms, as described in de¬ 
tail in Ref. [6] The remaining IR divergences are 


cancelled by the analogous singularities in the soft 
and collinear part of the real gluon emission cross 
section. 

The most difficult integrals arise from the IR- 
divergent pentagon diagrams with several massive 
particles. In Refs. [5,6] we calculated the pen¬ 
tagon scalar integrals as linear combinations of 
scalar box integrals using the method of Ref. [12]. 
For the gg initiated process we also used the 
method of Ref. [10] and found perfect agreement 
between the results of the two methods. The vir¬ 
tual corrections to the gg initiated process have 
an additional complication with respect to the qq 
case because of the presence of pentagon tensor 
integrals with rank higher than one. Pentagon 
tensor integrals can give rise to numerical insta¬ 
bilities due to the dependence on inverse pow¬ 
ers of the Gram determinant. The Gram deter¬ 
minant vanishes when two momenta become de¬ 
generate, i.e. at the boundaries of phase space. 
These are spurious divergences, which cause se¬ 
rious numerical difficulties. We have used two 
methods to overcome this problem and found mu¬ 
tual agreement within the phase space integration 
statistical uncertainty: 

• Impose kinematic cuts to avoid the phase 
space regions where the Gram determinant 
vanishes. Then apply an extrapolation pro¬ 
cedure from the numerically safe to the nu¬ 
merically unsafe region. 

• Eliminate all pentagon tensor integrals by 
cancelling terms in the numerator against 
the propagators wherever possible, after in¬ 
terfering the pentagon amplitude with the 
Born-matrix element. The resulting expres¬ 
sions are very large, but numerically stable. 

2.2. Real Corrections 

The real corrections are computed using the 
phase space slicing (PSS) method, in both the 
double [13] and single [14] cutoff approaches. In 
both approaches the IR region of the tth+g phase 
space where the emitted gluon cannot be resolved 
is defined as the region where the gluon kinematic 
invariants: 

Sig 2 Pi * Pg 2 htj Eg ( 1 {3i COS 0/g ) 


( 2 ) 
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become small. Here pi is the momentum of an 
external quark or gluon (with energy £)), Pi = 
'A — rrii/E'f, p g is the final state gluon momen¬ 
tum with energy E g , and 9i g is the angle between 
Pi and p g . In the IR region the cross section is 
calculated analytically and the resulting IR diver¬ 
gences, both soft and collinear, are cancelled, af¬ 
ter mass factorization, against the corresponding 
divergences from the 0(a s ) virtual corrections. 

The single cutoff PSS technique defines the IR 
region as that where 

Sig Smin , (P) 

for an arbitrarily small cutoff s m i n . The two cut¬ 
off PSS method introduces two arbitrary param¬ 
eters, 6 S and 5 C , to separately define the IR soft 
and IR collinear regions according to: 

rp 6 ,^ ci. 

Eg < —-— soft region , 

(1 — cos 6i g ) < S c collinear region . (4) 

In both methods, the real contribution to the 
NLO cross section is computed analytically below 
the cutoffs and numerically above the cutoffs, and 
the final result is independent of these arbitrary 
parameters. With this respect, it is crucial to 
study the behavior of cr NLO in a region where the 
cutoff(s) are small enough to justify the analytical 
calculations of the IR divergent contributions to 
the real cross section, but not so small as to cause 
numerical instabilities. 

This is the first application of the single cut¬ 
off phase space slicing approach to a cross section 
involving more than one massive particle in the 
final state. The numerical results of both meth¬ 
ods agree within the statistical errors. In Ref. [7], 
the dipole subtraction formalism has been used 
to extract the IR singularities of the real part. 
The agreement between these three very different 
treatments of the real IR singularities represents 
a powerful check of the corresponding NLO cal¬ 
culations. 

3. tth production at the Tevatron 

For pp collisions at y / s = 2 TeV, more than 95% 
of the tree-level cross section comes from the sub¬ 
process qq —> tth , while the gg and [q, q)g initial 
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Figure 1. Dependence of <Jlo,nlo(jpp —► tth) on 
the renormalization/factorization scale p, at the 
Tevatron ( v / i= 2 TeV), for M h = 120 GeV [5,6], 



M h (GeV) 

Figure 2. Dependence of ct lo ,nlo(pp —> tth) on 
Mh, at y/s = 2 TeV, for p = to * [5,6]. 

states are numerically irrelevant. Therefore, in 
[5,6], when calculating a NLO (pp —> tth ) of Eq. (1), 
we only included the qq — > tth channel, summed 
over all light quark flavors. The NLO inclusive 
total cross section is shown in Figs. 1 and 2, 
as a function of the renormalization/factorization 
scale p and as a function of the Higgs boson mass 
Mh, respectively. As expected, the NLO cross 
section at the Tevatron shows a drastic reduc¬ 
tion in the scale dependence from the lowest or¬ 
der prediction. The NLO corrections reduce the 
total cross section by a factor of 0.7 — 0.95 for 
rrit < p < 2mt- Only for p > 2mt + mh is the NLO 
cross section larger than the Born prediction. 
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Figure 3. cr NLO , NNLO for SM Higgs production 
processes at the Tevatron (y/s = 2 TeV) as a 
function of Mh- For pp —> tth, the renor¬ 
malization/factorization scale is varied between 
m t + M h /2< fj,<A(m t + M h / 2 ). 

To conclude, we want to summarize the state 
of the art of existing NLO and next-to-NLO 
(NNLO) calculations for the main SM Higgs bo¬ 
son production processes at the Tevatron. To 
this purpose, in Fig. 3 we show the renormaliza¬ 
tion/factorization scale dependence of the inclu¬ 
sive total cross sections at NLO (NNLO in case of 
gg —> h) for the main SM Higgs production pro¬ 
cesses. We only omit bbh associated production, 
since no complete NLO result is yet available. 
The renormalization/factorization scales are var¬ 
ied in the range 714/2< g<2Mh for all processes 
except pp — > tth , where we varied it in the range 
rn t + M h /2< g< 4 (m t + M h /2). The QCD NLO 
calculations for the qq —> Wh, Zh [15], qq —> qqh 
[16] and qq —> tth [5-7] processes provide reliable 
predictions for the inclusive total cross sections 
at the Tevatron. However, the NLO corrections 
to the gg —> h cross section [17] are large (up 
to ~ 100%), and the NNLO corrections, that re¬ 
cently became available [18,19], are crucial to ob¬ 
tain reliable predictions [20]. 

4. tth production at the LHC 

The associated production of tth at the LHC 
with yfs = 14 TeV is dominated at the parton level 
by gg —> tth , and in determining a NLO (pp —> tth), 
the 0(a s ) corrections to the sub-process gg —> 
tth provide the largest contribution. However, 
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Figure 4. Dependence of cr LO NLO (pp —» tth) on 
the renormalization/factorization scale g, at the 
LHC ( v / s = 14 TeV), for 714 = 120 GeV [9]. 



the other partonic channels cannot be neglected, 
and in calculating cr NLO {pp —> tth) of Eq. (1), 
we included all sub-processes, initiated by gg, qq, 
and (q,q)g [9]. 

The scale dependence of <J NLO (pp —► tth) at 
yfs = 14 TeV is also strongly reduced compared 
to the LO result, as can be seen in Fig. 4. The 
variation of <t nlo (pp —> tth) with the SM Higgs 
mass is shown in Fig. 5, for g = 2m t + Mh- At 
the LHC, the NLO QCD corrections enhance the 
LO cross section by a factor of 1.2—1.4 over 
the entire range of values of the renormaliza¬ 
tion/factorization scale shown in Fig. 4, and over 
the relevant Higgs boson mass interval shown in 
Fig. 5. 

Finally, as we did for the Tevatron in Fig. 3, we 
summarize in Fig. 6 the state of the art of existing 
NLO and next-to-NLO (NNLO) calculations for 
the main SM Higgs boson production processes 
at the LHC. 


5. Conclusion 

The inclusive cross section for pp, pp —> tth has 
been calculated at NLO of QCD, and the theoret¬ 
ical uncertainty on its prediction has been drasti¬ 
cally reduced to 10-15% for the Tevatron and to 
15-20% for the LHC. 

The availability of this NLO result as well as 
the recent calculation of the NNLO correction to 
inclusive Higgs production, allows us to summa- 
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M h (GeV) 

Figure 5. Dependence of <j lo ,nlo{pp —* tth) on 
Mh, at \fs= 14 TeV, for g = 2mj + Mh [9]. 



(GeV) 


Figure 6. a NLOtNNLO for SM Higgs produc¬ 
tion processes at the LHC (y/s = 14 TeV) as 
a function of Mh- For pp —> tth, the renor¬ 
malization/factorization scale is varied between 
m t + Mh/2<p<A(m t + Mh/ 2 ). 

rize the state of the art of higher order calculation 
for Higgs boson production cross sections at both 
the Tevatron and the LHC. The theoretical pre¬ 
dictions for Higgs boson production at the Teva¬ 
tron and the LHC are well understood at NLO, 
and except for the gg —> h process, for which the 
NNLO corrections are mandatory, the arbitrary 
scale dependence at NLO is small. 
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